A chiral crystal in cold QCD matter at intermediate densities? 
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The analogue of Overhauser (particle-hole) pairing in electronic systems (spin-density waves 
with non-zero total momentum Q) is analyzed in finite-density QCD for 3 colors and 2 flavors, 
and compared to the color-superconducting BCS ground state (particle-particle pairing, Q=0) . The 
calculations are based on efi'ective nonperturbative four-fermion interactions acting in both the 
scalar diquark as well as the scalar-isoscalar quark-hole (V) channel. Within the Nambu-Gorkov 
formalism we set up the coupled channel problem including multiple chiral density wave formation, 
and evaluate the resulting gaps and free energies. Employing medium-modified instanton-induced 
't Hooft interactions, as applicable around fiq ~ 0.4 GeV (or 4 times nuclear saturation density), we 
find the 'chiral crystal phase' to be competitive with the color superconductor. 



I. INTRODUCTION 
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^ , The understanding of QCD under extreme conditions is among the main frontiers in strong interaction physics. 
In particular, the finite-density and zero-temperature regime has reattracted considerable attention lately, after it 
has been realized that early perturbative estimates for color-superconducting gaps at large chemical potential are 
exceeded by up to two orders of magnitude towards smaller densities Such BCS-type pairing energies are in 

fact comparable to the ('constituent') quark mass gap in the QCD vacuum, Mq ~ 0.35-0.4 GeV, and have triggered 
new interest in observable consequences of quark matter formation within the core of neutron stars (unfortunately, in 
^ ' high energy heavy- ion collision large entropy production renders the access to this regime unlikely). 
— , The focus on the occurrence of various superconducting phases is motivated by the standard BCS instability of 
■ the Fermi surface for arbitrarily weak particle-particle (p-p) interactions. Under certain conditions, however, the 
' particle-hole (p-h) channel might also become competitive. Here, a kinematic singularity in the corresponding Greens 
function only develops in (effectively) 1-t-l dimensional systems, and at a total pair momentum of Q — 2pp (pp: 
Fermi momentum), known as Peierls instability In higher dimensions it can nevertheless be relevant provided 
. the interaction is strong enough. One variant of p-h instabilities are 'spin-density waves' as originally proposed by 
Overhauser Q for specific electronic materials (for a review see The analogue in the context of QCD, so called 
'chiral density waves', has first been discussed by Deryagin et al.\^. Using perturbative one-gluon exchange (OGE) 
at asymptotically high densities it was shown that the Overhauser-type pairing prevails over the BCS instability in the 
CIhI iVc — > oo limit {Nc. number of colors). This is due to the fact that the BCS bound states, being color non-singlet, are 
^ dynamically suppressed by 1/Nc as compared to the (colorless) Overhauser ones. More recently, Shuster and Son |^ 
revisited this mechanism for finite Nc and including Debye screening in the gluon propagator. As a result, the chiral 
density wave dominates only for a very large number of colors, Nc — 0(10"^). These findings have been confirmed in 
an analysis of coupled BCS / Overhauser equations using different arguments . One concludes that instabilities in 
the p-h channel are not relevant for real QCD at asymptotic densities. 
' The situation, however, can be very different if the interaction strength between the quarks is substantially increased 
(to be referred to below as the strong coupling regime). A well-known example is the Nambu-Jona Lasinio (NJL) 
description of chiral symmetry breaking in the QCD vacuum (associated with the constituent quark mass gap and the 
build-up of the chiral condensate), which requires a (minimal) critical coupling to occur. At finite density, the same 
(attractive) interaction is operative in the scalar-isoscalar channel. It's coupling strength is in fact augmented by 
a factor of {Nc — 1) over the (most attractive) scalar diquark channel. On the other hand, geometric factors act in its 
disfavor: unlike the BCS gap, which uniformly covers the entire Fermi surface, the chiral density wave appears in form 
of 'patches', their number depending on the symmetry of the presumed crystal. The purpose of the present paper is 
to study the interplay between Overhauser and BCS pairing, including different crystal structures, within the strong 
coupling regime. The focus is thus on quark matter at intermediate densities, i.e., large enough for the system to be 
in the quark phase, but small enough to support nonperturbative interactions. This should roughly correspond to 
chemical potentials in the range ~ 0.4-0.6 GeV, translating into baryon densities of 3.5-12 po (where po—0.16 im~^ 
denotes normal nuclear matter density) .0 



^The chiral crystal phase we are investigating is not to be confused with another crystal phase discussed in 1980's related 
to p-wave pion condensation llll] and later interacting skyrmions |13]. Those works have addressed nuclear matter at lower 
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The article is organized as follows. In Sect. || we start by introducing the Nambu-Gorkov type matrix propagator 
formalism that will subsequently be applied to obtain the gap equations for the coupled BCS-Overhauser problem; 
special attention is given to the single-quark spectra in the Overhauser ground state. In Sect. Ill we solve these 
equations using the aforementioned variants of nonperturbative interactions, i.e., somewhat schematic NJL-type as 
well as microscopic instanton-induced forces, for the slightly idealized case of two massless flavors and three colors. 
In Sect. [V we summarize and discuss the relevance of our results for real QCD. 



II. NAMBU-GORKOV FORMALISM AND COUPLED GAP EQUATIONS 



A. BCS Pairing 

A standard framework to address multiple instabilities in interacting many-body systems is provided by the Nambu- 
Gorkov formalism. Here, propagators are constructed as matrices combining all potential condensate channels via 
off-diagonal elements (see, e.g., ref. [p^), which automatically incorporates the interplay/coexistence of the various 
phases. 

For the familiar BCS case one adopts the following ansatz for the full propagator: 

r (k kA-n)-( ^^^■T'^I-t) (cfcTC-fci)\ _/G(fco,fc,A) F{ko,k,A) \ 
V (c-fei 4t) \(^-k[(^-ki) \F{ko,k,A) G{ko,-k,A)J 



The gap equation is then derived by formulating the pertinent Dyson equation, 

Gbcs 

which has the formal solution 



Go' A 



(2) 



Go ^Gq 1 - AA V -A Go 

where 



Gbcs - ZT ( ) ' (3) 



A = H)a,, I ^Fipo,p,A) (4) 



represents the (off-diagonal) 'selfenergy' contribution induced by p-p pairing (with an appropriate 4-fermion coupling 
constant app), and 

^0 = I 

ko - Ck + lOe^ 

Go = , ^ (6) 

are the free particle propagator and its conjugate at finite chemical potential (with = ojfc — fiq and infinitesimal 
(Jgj, = \S\ sgn(efc) according to the sign of e^). Inserting the expression for the anomalous Greens hmction from eq. (|^), 

[ko - ffc + «c>ej(fco + e/c + «c>ej - A^ 

into the definition of A, eq. (|^), yields the gap equation. Notice that the pole structure of F{k) always ensures a 
nonvanishing contour for the energy integration. 



densities, in which the chiral condensate < gg > is only slightly perturbed from its vacuum value and basically uniform in 
space, while the periodic structure is driven by pion fields. 
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B. Overhauser Pairing 



On the same footing one can analyze pairing in the particle-hole channel at finite total pair momentum Q. In the 
mean-field approximation (MFA)^, the full Greens function and Dyson equation in the presence of a single stationary 
wave take the form 



Govhiko,k,Q,a; fiq) = 



1 




f G{ko,k,Q,a) S{ko,k,Q,a) 
S{ko,k,Q,a) G{ka,k + Q,Q,a) 



Gn 



which has the formal solution 



S{ko,k,Q,a) 



(fco -ek+ i5ek){ka - (k+Q + i^e^+Q) - ^2 



and gives the ensuing gap equation from the definition of the pairing 'selfenergy', 



(2^ 



Notice that here the energy contour integration receives nonvanishing contributions only if 

^P+Q - cr^ < , 



(8) 



(9) 



(10) 



(11) 



which means that the two poles in po have to be in distinct (upper/lower) halfplanes, i.e., one particle (above the 
Fermi surface) and one hole (below the Fermi surface) are required to participate in the interaction. This condition 
reflects on the particle-hole symmetry caused by the nesting of the Fermi surface in the presence of the induced wave. 
We stress again that an important difference to the BCS gap equation resides in the fact that (for 2 or more spatial 
dimensions) one is not guaranteed a solution for arbitrarily small coupling constants since the p-h Greens function 
S does not develop a kinematic singularity (as mentioned in the introduction this is very reminiscent to the QCD 
vacuum case of particle- antiparticle pairing across the Dirac sea). 

At finite densities, the formation of a condensate carrying nonzero total momentum Q is associated with nontrivial 
spatial structures, i.e., crystals, characterized by a 'lattice spacing' a — 2t:/Q. In three dimensions a more complete 
description thus calls for the inclusion of additional wave vectors. In general, the p-h pairing gap can be written as 



(12) 



where the Qj correspond to the (finite) number of fundamental waves, and the summation over \n\ > 1 accounts for 
higher harmonics in the Fourier series. The matrix propagator formalism allows for the treatment of multiple waves 
through an expansion of the basis states according to 



/ (cfcT 4t) (cfct (cfcT •••\ 



G^ 



(Cfc+Qa 4t) (Cfc+Qa i^k+Q^l cl^c 

{Ck+Qyl 4|) {Ck+Qyl ci+Q^l) (Cfc+Q^ 4+C 



V 



(13) 



In practice the expansion has to be kept finite. The possibility of additional BCS pairing is straightforwardly incor- 
porated into eq. (O) by extending the latter with the off-diagonal states from eq. (n^). 



■^This approximation is equivalent to the weak coupling approximation in band structure calculations where higher intra-band 
mixing is suppressed. 
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In what follows we will consider up to riw = 6 waves in three orthogonal directions with = Qy = Qz and n = ±1, 
characterizing a cubic crystal through three standing waves with the fundamental modes (for simplicity we will also 
assume the magnitude of the various Overhauser condensates to be equal, i.e., Uj = <j). The important new features 
that arise through introducing additional states become already apparent in the simplest extension to 2 condensates. 
In this case one has for the (coupled) gap equation(s) 



d^p -<JxG^\p + Qy) 



(27r)4 Go i(p) G^\p + Qx) G^\p + Qy) - al G^\p ^ Qy) - al G^\p^ Qx) 



^(■\ f -(Ta; Go{p+Qx) 

[ *J aph J _ ^2 ^^^^^ _ ^2 ^^^^^ [ ) 

(and an equivalent one for ay by interchanging x <-* y). For finite Oy additional possibilities for the nonvanishing of 
the energy contour integration appear through an extra zero of the third-order polynomial in the denominator of the 
full propagator in eq. (IIj) . This enlarges the integration region and can be interpreted as interference effects between 
the patches (or waves). Diagrammatically this can be understood as an additional insertion of the Oy condensate on 
a particle (or hole) line of energy Cp. Note that a priori it is not clear whether such interferences are constructive or 
destructive, that is, give a positive or negative contribution to the right-hand-side (r/is) of eq. (|lj). 

In the propagators Gq contributions from antiparticles have been neglected. This should be a reasonable approx- 
imation in the quark matter phase at sufficiently large /ig, i.e., after the usual (non-oscillating particle-antiparticle) 
chiral condensate has vanished. At the same time, since we base our analysis on nonperturbative forces, the applicable 
densities are bounded from above. Taken together, we estimate the range of validity for our calculations to be roughly 
given by 0.4 GeV ^ /ig ^ 0.6 GeV. This coincides with the regime where, for the physical current strange quark mass 
of iris — 0.14 GeV, the two-flavor superconductor might prevail over the color-flavor locked (CFL) state so that our 
restriction to Nf = 2 is supported. 



C. Spectrum in the Overhauser Case 



The poles of the mean-field propagators discussed above provide the quasiparticle excitations in both the BCS and 
Overhauser case. In the former, the spectrum consists of gapped particles and holes. In the latter case, the physical 
interpretation is rendered more subtle by the presence of a standing wave. To keep the analysis transparent, we will 
discuss analytic results in 1+1 dimension and proceed to a numerical evaluation in 3-1-1 dimensions. 

In 1-1-1 dimension, the quasi-particle excitations following from the pole condition for the propagator in the Over- 
hauser case, eq. (^), have energies 

1 



e± = ■^{<^k + £k-Q) ± Y (cfc - Efe-g)^ + CT^ (15) 

with ek — \k\ ~ iiq and Q = 27r/a. This spectrum can be understood if we note that the quarks are moving in a 
self-induced potential V{x) = — 2o-cos (Qx) through the stationary wave. Indeed, in the presence of such a potential 
the spectrum is banded with |fc| < ir/a representing the first Brillouin Zone (BZ-1). In weak coupling the spectrum 
is mostly free except at A: = 0, dzn/a where band- mixing is large. For < Q we can ignore most of the band mixing 
except for the lowest one near the edge of BZ-1. Degenerate perturbation theory gives readily 



, (16) 



in agreement with (|5|). The 2x2 character of (|T^) follows from mixing between two bands. As foot-noted above, it is 
analogous to the MFA where the band-mixing is treated in the extended description commonly used in weak-coupling 
band-structure calculations. 
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^^=0.4GeV, Q=0.8GeV (4 waves), o=0 u =0.4GeV, Q=0.8GeV (4 waves), g=0.1 5GeV 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

k [GeV] k [GeV] 

FIG. 1. Dispersion relations of the various particle and hole branches using 4 waves in ±x and ±j/ direction with equal gaps 
and wave vector moduli Q = 2pF- The solutions are displayed for 3 angles of the momentum with respect to the a;-axis in the 
x-y plane. By symmetry, the solutions repeat itself within each octant (i.e., every 45°). The left panel shows the noninteracting 
case {a — 0), and the right panel has been obtained by setting a — 0.15 GeV. 

The quasiparticles of energy e_ are characterized by a standing wave ip-{x) w cos (nx/a), and those of energy e+ are 
characterized by a standing wave ip+i^) ~ sin(7ra:/a) near the edge of the Brillouin zone. The energy is substantially 
lowered by the standing wave ip-i^) with a probability density in opposite phase to the potential. The standing 
wave ijj+ix) corresponds to a probability density in phase with the potential, hence substantially more expensive 
energetically. At the edge of the zone, the two states are gapped by 2cr. Clearly, the lowest energy state is reached 
by filling only those states corresponding to that is by setting the Fermi energy at the gap. The ensuing state is 
an insulator. 

In higher dimensions, the band-mixing becomes more intricate. However, in the weak-coupling approximation and 
for Fermi momenta in the vicinity of Q/2, higher intra-band mixing is small and we may just use the extended band- 
structure description which is equivalent to our mean- field treatment. The quasiparticle spectra ej, {j = 0, . . . ,nm) 
follow numerically from the poles of the propagator.^ Fig. |^ shows an example of 4 waves in ±x and ±j/ directions 
for the canonical value of the wave vector, Q = 2pp^ with (a^ = > 0) ^-nd without {ax = Uy = 0) interactions. 
One clearly recognizes the formation of the gap close to the degeneracy point (level crossing) in the non-interacting 
case, which happens in the vicinity of the Fermi surface. For Q — pf (hg- ||) energy gain arises from an appreciable 
pushing down of the lowest level which lies rather deep within the Fermi sea. This is (partially) counteracted by an 
upward push of the upper branch which also corresponds to occupied states. Thus one expects the most beneficial 
configuration to be when the Fermi surface lies in between split levels. 



^Again, refers to the number of plane waves retained in a{x) with nm/2 being typically the number of standing waves. 
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^ =0.4GeV, Q=0.4GeV (4 waves), o=0 
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FIG. 2. Same as fig. ^ but for Q = pf. 

D. Energy Budget and Periodicity 

Solutions of the gap equations correspond to extrema (minima) in the energy density with respect to {wrt) the gap 
a. However, solutions may exist for several values of the wave vector Q. To determine the minimum in the latter 
quantity, one has to take recourse to the explicit form of the free energy density. In MFA, 



(17) 



where V3 is the 3-volume. The first contribution removes the double counting from the fermionic contribution in 
the mean- field treatment. Retaining only the particle contribution [i.e., neglecting antiparticles) , the free energy 
simplifies to 



E 



A 



BZ-l 



(18) 



with color/flavor coefficients C{Nc, Nf) which will depend on the concrete form of the pairing interaction. To avoid 
double counting the integration for the kinetic energy part is restricted to BZ-l as defined by the momentum regions 
[—Qj/2,Qj/2] (as well as \k\ < kp). This amounts to a folding of the various branches into BZ-l and enforces the 
explicit lattice periodicity onto the free energy. This point is illustrated in fig. ^ for the free case (cr = 0) along one 
spatial direction. For fixed chemical potential smaller wave vectors require the inclusion of an increasing number of 
branches to correctly saturate the available states within the Fermi sea (through a multiple folding until the Fermi 
surface is reached). E.g., for fjLq < 3Q/2 the lowest two harmonics with ± Qx suffice. Above, the next two higher 
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harmonics with ± 2Qx are necessary to encompass the occupied energy states within BZ-1 up to fig < 5(5/2, and 
so forth. When using additional waves in other spatial directions similar criteria hold (albeit more complicated due 
to nontrivial angular dependencies). 




FIG. 3. Free quasiparticle dispersion relations {aj = 0) for massless quarks in a crystal with periodicity in x-direction. 
The vertical dotted lines indicate the boundaries of the first {BZ-1) and second {BZ-2) Brillouin Zone (corresponding to 
< Qxfi and Qx/^ < |fca;| < Qx, respectively). Plotted are the 5 branches o;^ — \kx\, \kx ±Qx\ and \kx ±2Qx\ with the full 
lines characterizing their contributions to BZ-1. 

Solutions of the gap equations in general support different pairs of {cr, Q}; the combination that minimizes 
ri(/ig, Q, a) is the thermodynamically favored one. Note that the gap equations are not subject to explicit momentum 
restrictions since off-shell momenta of arbitrary magnitude can in principle contribute. 

As a simple example, ( p^ ) can be explicitly computed in 1+1 dimensions by recalling that the Fermi surface coincides 
with the gap, i.e., LVkp = fj.q — kp — Q/'2. Specifically, the contribution from the Fermi sea is 

nk^n=peF (^l-^ {Vl + e+eiry{^+^/lTe))'j , (19) 

with ^ = (j/Q and a density p — dfci^/vr where d is the overall degeneracy. In MFA, the induced standing wave is 
(j{x) = 2a cos{Qx), and the double counting in the Fermi sea is removed by 

The minimum value of a can be obtained in this case analytically by minimizing CI — Clkin + ^pot- 



III. NONPERTURBATIVE FORCES AND RESULTS 



For the actual calculations we need to specify the quantum numbers of the pairing channels. To do so we take 
guidance from low- (or zero-) density phenomenology encoded in effective 4-point interactions. In the particle-hole 
channel the strongest attraction is in the V channel given by (including exchange terms Q) 

^"mes = ^ , (21) 

whereas in the particle-particle channel it is believed to be the scalar diquark in the color-antitriplet channel (which, 
in fact, arises from a Fierz-transformation of (El])), 
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4., = g^2(jv^_i) (g^g75T2Alg) (g-r2A^75Cg-^) (22) 

(C: charge conjugation matrix, A^: antisymmetric color matrices, T2: SU(2)-flavor matrix). For practical use the 
effective vertices have to be supplemented with ultraviolet cutoffs. In the following we will consider two variants 
thereof and discuss the pertinent results for the coupled Overhauser/BCS equations. 



A. NJL Treatment 



In a widely used class of Nambu-Jona Lasinio models the ultraviolet behavior of the pointlike vertices is regulated 
by 3- or 4- momentum multipole formfactors (or even sharp 0-functions). We here employ a dipole form, 

— 2), for each in- and outgoing quark line with A = 0.6 GeV as a typical 'chiral' scale (variations within such 
parametrizations do not affect our qualitative conclusions in this section). The coupling constant A = 67 fm'^ is 
calibrated to a constituent quark mass of Mq = 0.4 GeV in vacuum. Note that there is no well-defined way of 
introducing density dependencies into the interaction. Since at finite fj,q the relevant quark interactions occur at the 
Fermi surface, a formfactor of type (|2^) implies the loss of interaction strength with increasing pp. 

This schematic treatment has been shown to yield robust results for 2-flavor BCS pairing with gaps A ~ 0.1 GeV 
at quark chemical potentials around 0.5 GeV [Q. Including now the p-h pairing as outlined in the previous section 
we find only rather fragile evidence for the emergence of chiral density waves (at fiq = 0.4 GeV): for wave vectors 
Qx < 0.150 GeV the rhs of the Overhauser gap equation supports solutions with gaps around ^ 5 MeV. The smallness 
of Qx in fact requires six waves (with k ± nQ^, n=l,2,3) to fill all states within the Fermi sphere. Somewhat more 
robust solutions are obtained when increasing the 4-fermion coupling constant. E.g., with a vacuum constituent 
quark mass of Mq — 0.5 GeV, which implies A = 73 fm^, the minimum solution emerges for ~ 0.2 GeV and 
cr ~ 20 MeV. However, the gain in the total free energy is very small: n = -1.2973 * lO'^ GeV"'' as compared 
to the free Fermi gas value of = —1.2969 * 10"'^ GeV"'' (to be contrasted with the BCS ground state for which 
^Bcsif^q = 0.4 GeV) ~ -1.375 * lO'^ GeV"" at a pairing gap of A ~ 0.13 GeV). 

We also checked that the incorporation of waves in other spatial directions does not lead to further energy gain. 



B. Instanton Approach at Finite Chemical Potential 



A more microscopic origin of effective 4-fermion interactions is provided within the instanton framework. In the 
finite-density context it has previously been employed to study the competition between the chiral condensate and 
two-flavor superconducting quark matter in refs. |14 ^. Let us briefly recall some elements of the approach. The 
starting point is the QCD partition function in instanton approximation. 



1 



n 

1=1 



dflj n{pi) e 



[det{i p - ifiqj4)]'^f , 



(24) 



where flj — {Qj, pj, zj} denote the collective coordinates (color, size and position) of the instanton solutions and 
n{pj) their individual weight. To extract effective quark interactions one reintroduces quark flelds in a way that is 
compatible with the fermionic determinant of the previous equation p5| , 



Ztnst Uj-q) = J Vil^Vip^ exp j d'^xip'' (i ^ - ipq"fi)ll^ 



dX± 
2n 



exp 



A±y± + iV± In 



N± 



1 



(25) 



and an additional auxiliary integration over A± has been introduced to exponentiate the effective (2A^^)-fermion 
vertices Y± . For two flavors the latter are given by 



A± Y± = A± 



d'^kid'^k2 d'^p id'^p2 



(27r) 



16 



{2T,)H^^\k^+p^-~k2-p2) 



X [^P^^\pl,-^^q)l±Ta •?^(fcl,Mg)^] [i^'^^\p2,~p.q)l±Ta T{k2,Pq) 



(26) 
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with 7-1- = (1 ± 75), flavor matrices = (f, i) and the instanton formfactors ^(p, Hq) = {p + inq)~ ^p{p, /iq)"*", which 
are matrices in Dirac space, adopting the notation of ref. [Q, i.e., = x^^a^ with cr^ = (±(t, 1). Since the 
fermionic determinant has been approximated by its zero- mode part, the formfactors are entirely determined by the 
Fourier-transformed quark zero- mode wave functions, 

4>i^a{p, ^iq) = j d^x e"*^''^ 4>I^A{x,flq) 

= V{P,f^q)^ XR.L (27) 

with (j){x,pLq) satisfying the Dirac equation in the background of an (anti-) instanton, 

pi^A ^lA^) = • (28) 

The explicit form of Lp^{p,^q) can be found in refs. p^,p[. 

As before we preselect the potential condensation channels as the scalar-isoscalar p-h and p-p ones which, after 
solving the (matrix) Dyson equation, yields the coupled gap equations 

A = j^gi^^ j 1^ Bip, Fip; A, a„Q,) (29) 

{ax + 5Sax)= — J j^-:^ A{p,Qx;fj.q) S^ip; A, aj , Qj) (30) 



{as,. - 10 Sa,) ^ y ^ A(p,g,;^,) 5,(p;A = 0,(T3,„Q,) . (31) 



Here, 5a j — aj — a^j denotes the difference in the Overhauser gaps for quarks of color 1,2 (cj) or color 3 (0-3 j ) which, 
respectively, do (eq. (pO|)) or do not (eq. (^) participate in the diquark pairing (once A 7^ 0) jl^,^. In principle, the 
wave vectors in the color-1,2 and color-3 channels could also be different when minimizing the free energy. However, 
the actual solutions for aj and CT3j turn out to be very close to each other even in the presence of large ECS gaps 
A so that the generically very smooth dependence on the Qj should not cause appreciable deviations between the 
two color sectors. Under our simplifying assumption that the momentum moduli \Qj\ (as well as the associated gap 
parameters aj) of the Overhaus er p airing are of equal magnitude the additional {2nio — 2) gap equations for the other 
p-h channels are equivalent to ( ^0^ and (^). The explicit form of the propagators is given by 

Sxip,A,aj,Qj) = -~A{p,Qx;^iq) (Jxip^Qx) Go{p + Qx) D{p,A,aj,Qj) (32) 
Fip,A,aj,Qj) = -Bip;,ig) A(p) Go(p) D{p,A,a,,Qj) (33) 

with 

D{p, A, a,,Q,) = [A{pfGo{p) - G^\p) + a,{p, Q, f G„{p + Q,)] . (34) 

j 

The functions A and B represent the (square of the) instanton form factors (normalized to one in vacuum) acting 
on each fermion line entering/exiting a vertex, and we have introduced the notation aj{p] Qj) = ajA{p; Qj), A(p) = 
AB{p). The integration variable A = A± plays the role of an effective coupling constant, which, however, is not 
a priori fixed. Rather, its value is found from minimization of the free energy via a saddle point condition, which 
reads Q 

N 

- = A(y+ + y_), 

j 

where (1+ + Y^) denotes the ground state expectation value of the interaction vertices with potential condensates. 
Thus the magnitude of the gaps itself governs the effective coupling to the instantons. In the present treatment the 
instanton density N/V is assumed to be constant The final result for the free energy at the minimum then becomes 



*This assumption is motivated by the observation that the free energy associated with the instanton vacuum (background) 
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InZ 



N 

= ri/cm(Mg) + y 111 



(N/V) 



(36) 



indicating that the potential (second) term favors small values for A, whereas the kinetic (first) term exhibits the 
usual decrease with increasing values for the gaps/condensates (and thus for A). We should also point out that flinq) 
is only determined up to an overall constant which is associated with the nonperturbative vacuum energy of about 
—0.5 GeV/fm'^ (or, equivalently, bag pressure P>0). This term is encoded in the scale dependence of the argument 
in the logarithm (note that A and N/V have different dimensions), which we have not assessed here since it is not 
relevant for our analysis (this will be the origin of positive values for ^l{nq) encountered below). 

Before we come to the numerical solutions of the gap equations let us recall the specific density dependence of the 
instanton formfactors, cf. fig. ^ At fixed energy (upper panel) the strength of the interaction is clearly concentrated 
at the Fermi surface: the falloff with three- mometum sets in only above p c:i pp. On the other hand, as a function 
of (Euclidean) energy the strength is reduced starting from p4 = (see lower panel). These features reflect that the 
instanton zero- modes, which mediate the interaction, operate across the Fermi surface, i.e., at zero energy but at 
3-momenta equal to the Fermi momentum. This behavior is quite distinct from the schematic (density-independent) 
NJL forces employed in the previous section. Already at this point one can anticipate the Overhauser pairing to be 
more competitive than in the NJL treatment. 




FIG. 4. In-medium instanton form factors as a function of three-momentum. Upper panel: at fixed energy p4 = for 
chemical potetnials jiq — 0, 0.3 GeV, 0.6 GeV; lower panel: for fixed /i, — 0.3 GeV and various energies. 

For the evaluation of the kinetic part of the free energy as given in eq. ( |l8| ) a complication arises from the fact that the 
the instanton formfactors are defined in Euclidean space. We therefore approximated the gaps entering into the integral 
for flkin by their zero-energy values retaining the 3-momentum dependence, i.e., A{p) ~ A(p) = AB{pi = 0,p), and 
equivalently for a. It turns out that this approximation is consistent in the sense that the resulting extrema in 



is rather large compared to interaction corrections arising in the finite-density quark sector. It is corroborated by explicit 
calculations in the 'cocktail model' of ref. , where the grand potential has been minimized explicitly over both the properties 
of the instanton ensemble as well as quark Fermi sphere: the resulting variations in the total N/V where found to be at the 
few percent level. However, the assertion of constant N/V can imply unphysical behavior when the system is driven towards 
very large or very small pairing gaps. 
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are in line with the solutions of the gap equations (which are solved in Euclidean space with the full 4-monientuni 
dependence of the formfactors) . Other choices for fixing p4 do not comply with this criterium. 

If not otherwise stated, the subsequent calculations have been performed for a total instanton density of N/V — 
1 fm"** which in vacuum translates into a constituent quark mass of Mq — 0.34 GeV. The search for simultaneous 
solutions to the gap equations (p9|), ( [3C| ) and (^ij) together with the selfconsistency condition on the coupling, cq. (|3^), 
is illustrated in fig. |5. The upper two curves represent the values for a and (73 that simulatneuosly solve the two 
Overhauser gap eqs. (30|) and ( |3l| ) at given BCS gap A (plotted on the abscissa). The lower curve indicates the values 
for a that solve the BCS gap eq. ( p9| ) for given A (using the value for as from the upper curve to fix the coupling A). 
Thus a coexistence state of Overhauser and BCS pairing would be signalled by the crossing of the two lower lines. 
As mentioned above no such state is found; the only physical solutions correspond to the crossing points of the upper 
(full or dashed) curve with the y-axis (Overhauser state with A = 0, a — finite) and of the lower (dashed-dotted) 
curve with the a;-axis (BCS state with a = (73 = and A = 0.225 GeV; the discrepancy of about 15% (less at higher 
/ig) with refs. |14[1] reflects the accuracy when neglecting antiparticle states). 




FIG. 5. Solutions of the coupled Overhauser and BCS equations with 2 waves at Qx ~ ±0.5 GeV. The full and dashed lines 
indicate, respectively, the values for aj:{A) and 0-3, ^ (A) which simultaneously solve the two Overhauser gap equations ( po[ ) and 
for given A. The dashed-dotted line corresponds to the points ax, A that solve the BCS gap equation (|2^). A nontrivial 
simultaneous solution to all 3 gap equations would be signalled by a crossing of the full and dashed-dotted curves. The two 
independent pure BCS and Overhauser solutions are marked by the squares on the x- and y-axis, respectively. 

The question then is which one is thermodynamically favored. The BCS solution is unique (A = 0.225 GeV) and 
gives a total free energy of ^Bcsil^q — 0.4 GeV) = 2.3 * 10~^ GeV^ (up to a constant which is not relevant here, as 
discussed above). 

The situation is more involved for the Overhauser configurations. Let us start with the 'canonical' case where the 
momentum vector of the chiral density waves is fixed at twice the Fermi momentum, Q = 2pF- In fig. |^ the resulting 
minimized free energy (corresponding to solutions of eq. (p0|)) is displayed as function of the number of included 
waves. The Overhauser solutions are not far above the BCS groundstate, with a slight energy gain for an increased 
number of waves. 

However, one can further economize the energy of the Overhauser state by exploiting the freedom associated with 
the wave vector Q (or, equivalently, the periodicity of the lattice). For Q > 2pp the free energy rapidly increases. On 
the other hand, for Q < 2pp more favorable configurations are found. To correctly assess them one has to include 
the waves in pairs \k ± Qj \ of standing waves (i.e., — 2, 4, 6, . . .) to ensure that the occupied states in the Fermi 
sea are saturated within the first Brillouin Zone (cf. fig. ||). The lowest- lying state we could find at Hq = 0.4 GeV 
occurs for one standing wave with Qmin — 0.5 GeV and a ~ 0.21 GeV with a free energy ~ 2.3 * 10~^ GeV''', 
practically degenerate with the BCS solution. In solid state physics the breaking of translational invariance in one 
spatial direction is typically associated with 'liquid crystals', i.e., 2-dimensional layers of (uniform) liquid separated 
by periodic spacings of a = 2-k/Q. In our case, a ^ 2.5 fm. The minimum in the wave vector is in fact rather shallow, 
as seen from the explicit momentum dependence of the free energy displayed in fig. ^ 
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FIG. 6. Dependence of the free energy (upper full line) and associated p-h pairing gap (dashed-dotted) on the number of 
included waves ('patches') with fixed 3-momentum modulus \Qj\ — 0.8 GeV for solutions of eqs. ( |30| ) and (|35|); the full line 
marks the value of the BCS ground state free energy that solves eqs. (^^ and (|35|). The results are for fiq — 0.4 GeV and 
N/V = 1 fm"*. 
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FIG. 7. Left panel: wave- vector dependence of the Overhauser free energy for one standing wave (fuU hne: Q.ZT , short-dashed 
line: ^^i^) in comparison to the BCS solution (long-dashed line: Q,f^^ , dotted line: 57^,^'^) at = 0.4GeV (we note again 
that the absolute values of the total free energies are only determined up to an overall (negative) constant (related to the 'bag 
presssure') which drops out in the relative comparison of the solutions); right panel: wave- vector dependence of the Overhauser 
pairing gap (full line) compared to the BCS gap (long-dashed line). 



Finally we confront in fig. g the density dependence of the free energies (upper panel) and pairing gaps (lower 
panel) in the minimum of the BCS and n^=2,6 Overhauser states. Again we see that over the applicable /x^-range the 
solutions are close in energy, with (almost) degenerate minima for the n^=2-'liquid crystal' and the BCS ground state 
at the lower densities of ps — 4po- Towards higher densities, where the gaps and thus the strength of the effective 
instanton interactions decrease, the BCS solution becomes relatively more favorable. This confirms once more that 
in the 3-f 1 dimensions the Overhauser pairing can only compete for sufficiently strong coupling. We should also note 
that for both the 71^=2- and riuj=6-case as displayed the optimal momentum vector of the pertinent standing waves 
stays at approximate values of Qj ~ 5pi?/4 and bpF/i, respectively. 
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FIG. 8. Chemical potential dependence of the total free energies (upper panel) and pairing gaps (lower panel) for the BCS 
(long-dashed lines) and Overhauser (solid and dashed-dotted lines) solutions. 



IV. CONCLUDING REMARKS 



Employing a standard Nambu-Gorkov (matrix) propagator approach we have performed an analysis of competing 
instabtilities in the particle-particle and particle-hole channel for three-color, two-flavor QCD at moderate quark 
densities. As an essential ingredient we used nonperturbative forces (strong coupling) and preselected the potential 
condensation channels with guidance from low-energy hadron phenomenology, i.e., both the diquark as well as the 
quark-hole pairing were evaluated in their scalar-isoscalar channels. The corresponding coupled gap equations do not 
seem do support simultaneous nontrivial solutions. Needless to say that our calculations might be further improved 
by including higher standing waves through a larger class of crystalline symmetries (polyhedron structures). Also, 
our mean-field approximation may be extended to account for higher intra-band mixings when > Q. 

This not withstanding, an important outcome of our calculation is that the individual (separate) solutions for the 
BCS and Overhauser ground states are quite close in energy, indicating the importance of particle-hole instabilities 
when using interaction strengths as typical for low-energy hadronic binding. This becomes particularly evident when 
comparing schematic NJL interactions with microscopic instanton forces: whereas the former are reduced in magnitude 
with increasing density (entailing a relative suppression of the Overhauser pairing), the density-dependent instanton 
form factors essentially preserve their strength at the Fermi surface rendering the Overhauser pairing competitive. 
Indeed, with a somewhat larger instanton-density of N/V = 1.4 fm~*, corresponding to a vacuum constituent quark 
mass of Mq — 0.41 GeV, the Overhauser solution reaches below the BCS one. On the other hand, if we were to 
minimally account for strange quarks, the instanton interaction (being a 6-quark vertex) in the ud sector would lose 
about 60% of its strength due to the reduced strange-quark mass {iris — 0.15 GeV as opposed to Mg ~ 0.45 GeV in 
vacuum) in closing-off the strange quark line. Some of this loss might be recovered once strange quarks themselves 
participate in the Overhauser pairing. Our findings are to be contrasted with earlier calculations based on perturbative 
OGE at high densities where an extremely large number of colors was required for the Overhauser pairing (i.e., a 
chiral density wave) to overcome the BCS instability. 

Finally, a remark about the relevance of our results for neutron stars is in order. Here, quark matter is believed to 
reside mostly in a mixed phase, with significant charge separation between quark and nuclear components. The quark 
core, however, also exhibits charge asymmetry due to the finite strange quark mass. Consequently, quark matter 
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should have appreciably different up- and down-quark chemical potentials, which suppresses the flavor-singlet diquark 
(ud) pairing On the other hand, such a suppression does not apply to the (flavor-singlet) particle-hole channels of 
type uu^^, dd~^, which suggests that isopin asymmetric quark matter provides additional favor to the chiral crystal. 
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